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SECTION  1 
INTRODUCTION 


This  report  considers  data  and  mathematical  models  for 
geological  materials  subjected  to  shock  loading.  Data  from 
laboratory  experiments  are  summarized  in  terms  of  pressure  (P) 
and  density  (p)  at  equilibrium  states.  An  ensemble  of  such  states 
over  a  range  of  pressures  constitutes  the  principal  Hugoniot  for  a 
material.  The  specific  internal  energy  (e)  at  these  states  can 
be  found  by  applying  the  Rank! ne-Hugoniot  equations  for  conserva¬ 
tion  of  mass,  momentum,  and  energy.  Data  are  also  given  for  states 
along  adiabatic  pressure  release  paths.  The  model  must  be  used 
cautiously  below  50  kb,  because  variations  in  density,  porosity, 
and  moisture  within  the  same  basic  material  from  different  loca¬ 
tions  can  dramatically  affect  the  low-pressure  properties.  The 
model  accounts  for  these  variations  if  the  empirical  constants  are 
properly  modified.  At  high  levels  of  pressure  and  specific  internal 
energy,  the  model  represents  a  perfect  gas  (gamma  law  gas,  i.e. 

P  =  (y-1)  pe)  with  {  =  1.5-  Under  conditions  of  very  low  density 
and  specific  internal  energy  above  the  vaporization  energy  (em) , 
the  model  represents  a  perfect  gas  with  y  =  1.1. 


Analytic  equations  of  state  are  provided  which  match 
approximately  the  main  features  of  the  data  over  the  range  men¬ 
tioned.  The  equation  of  state  aims  primarily  to  calculate  P, 
once  p  and  e  are  specified.  The  basic  form  of  the  high- 
pressure  model  is  due  to  Tillotson  (Reference  1-1).  However, 
deviatoric  stress/strain  relations  are  accounted  for  in  order  to 
provide  a  more  realistic  model  in  the  low-pressure  region.  The 
deviatoric  stress/strain  relations  are  prescribed  by  means  of  a 
variable-moduli,  perfectly  plastic  nodel ,  wi th  associated  plastic 
potential  flow  rule.  This  model,  when  used  wi th  a  Coulomb-type 
yield  criterion  depending  on  mean  stress,  produces  di latency  or 
plastic  volumetric  expansion.  This  tendency  opposes  the  plastic 
volumetric  compaction,  or  hysteresis,  which  enters  through  the 
variable  bulk  modulus.  Although  rigorous  proofs  of  uniqueness  are 
not  generally  available  for  such  an  analytic  model  (Reference  1-2), 
which  in  any  event  apply  only  in  the  low-pressure  range,  no  diffi¬ 
culties  in  obtaining  unique  solutions  have  been  encountered  in 
practice. 


The  major  modifications  to  previous  work  done  in  this 
area  (References  1-3  and  1-4)  include  accounting  for  (1)  the 
specific  energy  lost  or  gained  during  phase  changes,  (2)  changes 
in  bulk  and  shear  moduli  and  in  shear  strength  due  to  phase 
changes,  and  (3)  volumetric  compaction  due  to  irreversible  closure 


1 


I 


a 


of  cracks  or  pore  spaces.  Also,  specific  parameters  of  the  modal 
for  ?ach  material  are  selected  after  considering  both  Hugoniot  and 
release  i sen  trope  data. 

The  materials  considered  in  this  study  include  the 

following: 


NTS  Granite  (granodiori te) 

Cedar  City  Tonal ite 
Laramie  Anorthosite 
Banded  Mountain  Limestone 
Coconino  Sandstone 

Mountain  Home  Basalt  (dense  and  porous) 

Tuff 

Sal  t 


The  data  on  which  each  model  is  derived  are  obtained 
primarily  from  laboratory  experiments  on  samples  whose  maximum 
dimension  is  1  to  2  in.  Site  surveys  are  also  considered  for 
general  rock  classification,  in  situ  bulk  density,  porosity, 
water  content,  seismic  velocity,  crack  patterns,  and  geologic 
layering.  It  is  common  in  such  a  survey  to  find  variations  of 
as  much  as  10  to  15  percent  in  as  basic  a  property  as  the  density. 
Other  properties,  such  as  porosity  or  seismic  velocity,  may  easily 
vary  more  than  a  factor  of  two  within  a  geological  area  much 
smaller  than  the  region  considered  in  a  ground  motion  calculation. 
Judgment  is  a  necessary  ingredient,  therefore,  in  selecting  repre¬ 
sentative  values  for  each  of  the  parameters  measured  in  these 
surveys. 


Because  the  tests  included  in  a  site  survey  are  performed 
at  low  pressure,  one  must  turn  to  laboratory  results  to  obtain  the 
high-pressure  dynamic  properties  of  the  rock.  However,  the  site 
surveys  are  useful  for  identifying  the  specimens  which  are  most 
representative  of  the  in  situ  medium  and,  based  on  the  homogeneity 
of  the  site,  for  determining  how  extensive  the  laboratory  programs 
must  be  to  describe  the  rock  in  the  field.  At  mean  stress  levels 
less  than  a  few  kilobars,  the  triaxial  and  hydrostatic  tests  pro¬ 
vide  most  of  the  data  necessary  to  determine  the  values  of  the 
parameters  in  the  constitutive  model.  The  hydrostat  provides  the 
relation  between  mean  stress  and  relative  volume,  and  triaxial 
tests  are  used  to  determine  the  shear  modulus  and  yield  criterion. 
Loading  and  unloading  cycles  in  both  types  of  tests  are  necessary 
to  measure  the  extent  of  hysteretic  compaction  in  porous  media  at 
various  loading  stresses.  Wave  speed  measurements  and  low-pressure 
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Hugoniot  data  provide  a  useful  check  on  these  properties.  At 
high  pressures,  the  Hugoniot  data  provide  a  measure  of  the 
material ‘s  elastic  limit  and  its  behavior  up  to  energy  densities 
at  which  the  Thnmas-Fermi  theory  (i.e.,  perfect  gas  with  y  =  1.5) 
can  be  applied,  release  Isentropes  from  mul tiki  lobar  shock 
states  help  complete  the  description  of  hysteretic  compaction. 
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SECTION  2 

GENERAL  EQUATIONS  OF  STATE 


The  general  equations  of  state  are  divided  into  two 
parts.  The  first  part  describes  calculation  of  the  pressure  or 
mean  stress  and  the  second  part  describes  calculation  of  devia- 
toric  stresses  in  the  solid  phase.  In  the  fluid  state,  the 
deviatoric  stresses  vanish  and  only  the  first  part  of  the  EOS  is 
necessary. 


CALCULATION  OF  PRESSURE  OR  MEAN  STRESS 


Loading 

The  mean  stress  or  pressure  is  expressed  as  the  sum  of 
two  independent  parts:  the  fluid  pressure  (Pf)  and  the  solid 
pressure  (Ps) . 


P  -  Pf+Ps 


(2-1) 


Since  the  pressure  in  a  fluid  is  assumed  to  be  independent  of  its 
previous  history,  Pf  is  a  function  of  only  specific  internal 
energy  (e)  and  density  (p).  The  equations  for  calculating  the 
fluid  pressure  are  a  modification  of  the  model  first  proposed  by 
J.  Tillotson  (Reference  2-1),  for  metals  subjected  to  strong 
shocks: 


Pf  “ 


exp  Z  pe* 


(2-2) 


where 


a2’  k’  eo  =  ^mP'r'ca^  constants 

=  (p  is  reference  density) 

pQ 
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Cr i terion 


Material  State 


u  >  0 


Liquid  or  compressed  gas 


u  <  0 


Fxpanded  gas 


0 

a;':  = 

e 


e  <  e  Cold  sol  id 

m 

e  >  e  Fluid  with  specific  internal 

energy  greater  than  that 
required  to  begin  the  phase 
change  (em)  from  the  solid  state 


The  specific  internal  energy  at  which  the  phase  change  is  initiated 
depends  on  the  compression  n 


e  = 
m 


e  (1  +  fn)  <  e 
mo  mm 


where 

e^  »  Specific  internal  energy  required  to  initiate  a 
phase  change  at  normal  density 

e  =  Specific  internal  energy  at  which  a  phase  change 
11  is  completed 

f  =  Empirical  parameter 

The  role  of  exp  z  is  to  provide  a  smooth  transition  from  the  com¬ 
pressed  state  to  that  of  an  ideal,  low  density  gas.  The  reason 
for  defining  a  new  variable,  e*,  to  represent  the  internal  energy 
is  to  account  for  the  energy  absorbed  or  released  during  a  phase 
change  which  does  not  affect  pressure  or  density.  Thus,  on  load¬ 
ing,  there  is  no  contribution  to  P  from  Pf  until  e  =  em.  As 
loading  continues  and  e  >  e  ,  it  is  shown  below  that  the  contribu¬ 
tion  of  Ps  rapidly  decreases  and  is  replaced  by  Pf.  Conversely, 
as  unloading  proceeds  from  a  fluid  state,  the  main  contribution  to 
P  is  from  Pf  until  e  =  em.  If  unloading  continues  such  that 

e  <  e„.  Pc  becomes  the  dominant  contributor, 
m’  a 

Under  conditions  of  low  density  and  high  internal  energy, 
Equation  2-2  effectively  reduces  to 


Pf  =  aj  pe* 


(2-3) 


5 


R-71 34-2283 

This  is  the  familiar  equation  for  a  low  density  perfect  gas  in 
which  aj  =y  ~  l  where  y  is  the  ratio  of  the  specific  heats  at 
constant  pressure  and  volume.  For  purposes  of  the  calculations, 
it  is  assumed  that  a^  is  a  constant,  even  though  y  is  known 
to  vary  (Reference  2-<.)  with  P,  0,  and  -{.  To  minimize  the 
error,  an  average  value  of  aj  ir,  the  ranje  of  interest  is 
selected.  The  values  a£  and  b  are  choset  to  match  experimen¬ 
tal  Hugoniot  data  at  pressures  of  t'le  order  of  500  kb,  while 
e^  Is  selected  so  that  the  model  approximates  the  Thomas-Fermi- 
Dirac  description  of  highly  compressed  material  at  high-energy 
densities. 

Calculation  of  the  mean  stress  in  the  solid,  Ps  in 
Equation  2-1,  is  based  on  J.  B.  Walsh's  concept  (References  2-3  and 
2-4)  that  the  effective  values  of  the  elastic  parameters  differ 
from  the  intrinsic  values  for  the  consolidated  material  due  to  the 
presence  of  cracks  and  pores.  The  hydrostatic  loading  behavior  of 
a  material  containing  cracks  and  pores  is  then  described  by  a  bulk 
modulus  which  varies  in  the  following  way: 

3P  /  v 

T-5-  =  K  -  (K  -  K  )  exp  (2-4) 

3u  m  '  m  o  K  \\i-  J 

where 


u 

u* 

Integrating  Equation  2-4  between  the  limits  0  and  p  +  Be  (B 
is  the  coefficient  of  volumetric  thermal  expansion  per  unit 
energy),  leads  to  the  following  equation  for  Pg: 

Ps  =  Km  <„  .  Be)  -  ;kb  -  Ko)  „*  [l  -  exp  (2-5) 

Initially,  p  +  Be  =  0  and  the  bulk  modulus  is  K  ,  the  bulk 
modulus  of  the  rock  matrix  with  included  voids.  As  p  increases 
the  voids  close,  and  the  bulk  modulus  approaches  the  intrinsic 
value  Km.  The  function  of  the  empirical  constant  p*  is  to 
control  the  rate  at  which  the  solid  bulk  modulus  asymptotically 


=  Initial  bulk  modulus 
=  Intrinsic  bulk  modulus 


=  Elastic  component  of  excess  compression  = 
=  Empirical  constant 
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approaches  its  maximum  value.  In  selecting  a  value  of  .j*  for  a 
specific  material,  it  is  helpful  to  remember  that  jPs/ will 
approach  Kn]  faster  as  t.*  is  made  smaller.  This  formulation 
assumes  that  the  mean  stress  depends  only  on  the  elastic  component 
of  volume  change.  Inelastic  volume  change  due  to  plasticity  and 
strain  rate  effects  are  specifically  excluded  from  computation  of 
pressure. 

Ps  depends  on  e  in  two  ways.  First,  u  is  augmented 
by  the  thermal  expansion,  Be,  so  that  a  hot  solid  is  at  a  higher 
pressure  than  a  cold  one  at  the  same  density.  Also,  since  the 
intrinsic  bulk  modulus  of  most  rocks  decreases  with  increasing 
temperature  or  energy  (Reference  2-5),  Km  in  Equation  2-5  is 
replaced  by 

Km  =  Kmax['  '  (C  r)\ 
where 

K  =  Intrinsic  bulk  modulus  at  room  temperature 
max 

c  =  Empirical  constant 

Thus,  the  solid  contribution  to  the  i.'fective  bulk  modulus  becomes 
zero  when  the  phase  change  is  initiated.  The  effect  of  causing 
the  bulk  modulus  to  decrease  with  increasing  specific  internal 
energy  is  to  make  a  portion  of  the  principal  Hugoniot  to  be  con¬ 
cave  to  the  y-axis,  as  is  illustrated  in  Figure  2-1.  In  selecting 
values  of  the  parameters  c  and  f  for  a  specific  material,  it  is 

helpful  to  remember  that  point  m  moves  toward  the  origin  if  c 

is  increased,  if  e  is  decreased,  or  both, 
mo 

The  model  described  above  is  adequate  for  materials  whose 

initial  porosity  is  less  than  about  5  percent.  For  materials  with 

greater  initial  porosity,  the  very  low-pressure,  solid-phase  hydro- 
stat  is  modified  as  shown  in  Figure  2-2(a).  The  material  is 
assigned  a  reference  density  Pref»  which  is  assumed  to  be  the 
density  at  zero  pressure  along  a  master  unloading  curve  as  shown 
in  Figure  2-2.  A  new  variable  u  is  defined,  which  enables  Uiu 
P/u  relation  on  loading  to  be  calculated  from  the  master  unloading 
relation.  Thus 


-  a  •  e  ■  -nVfnT  *  0 

i  mo  J 


A 


where 
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u  <  p 


U1  <  M  <  U2 


u2  <  v 


n  *  i>  -  m0 

p  =  +  u2j  exp 

u  *  P 


(u2  -  u3) 


The  empirical  parameters  pj  ,  p2,  and  p^  have  the 

following  physical  interpretation: 

p.  =  Excess  compression,  at  which  the  material  matrix 
begins  to  break  down 

u2  =  Excess  compression,  at  which  the  virgin  loading 
curve  joins  the  master  unloading  curve 

u,  =  p,  -  u 


Unloading 

In  the  present  model,  primary  consideration  during 
unloading  is  given  to  hysteresis,  the  irreversible  compaction 
which  may  occur  during  a  cycle  of  hydrostatic  loading  and  unload¬ 
ing.  As  an  example,  experimental  data  (Reference  2-6)  indicate 
that  UTS  granite  compacts  irreversibly  to  about  1.005  times  its 
initial  density  when  it  is  subjected  to  pressures  between  5  and 
50  kb  and  then  unloaded.  If  the  peak  loading  pressure  is  less 
than  5  kb,  the  hysteretic  compaction  is  less  than  0.005,  but 
because  the  unloading  data  are  limited,  the  functional  relation¬ 
ship  between  the  permanent  compaction  and  peak  pressure  in  this 
regime  is  based  on  hypothesis. 

The  mathematical  procedure  for  representing  hysteresis 
is  an  extension  of  the  variable  modulus  method  used  for  the  load¬ 
ing  relationship.  If  p  is  greater  than  or  equal  to  the  maximum 
excess  compression  ( p  )  previously  experienced  by  the  material, 

Ps  is  calculated  by  Equation  2-5,  which  is  referred  to  as  the 
virgin  loading  curve.  However,  if  u  <  pm,  the  material  is 
on  an  unloading/ reloading  path.  The  permanent  compaction  or 
"set"  resulting  from  loading  to  pm  is  then  calculated.  This 


set  represents  a  new  point  along  the  p-axis,  namely 


at  which 


10 


.sJfray’1*  ».  -  ..  _ 


R-71 3^-2283 

the  mean  stress  returns  to  zero.  The  unloading  path  has  the 
same  value  of  pressure,  Pm,  as  the  virgin  loading  curve 
when  u  =  Thus,  If  Kg,  and  u*  are  assumed  to  be  constant. 
Equation  2-6  can  be  inverted  to  obtain  a  new  value  of  Kq, 
namely 


K‘  =  K  P  -  k(u  -ii  +  Be)!  x 

o  m  I  m  in'  m  z  / 1 

M1  •  exp  [-  K  -  uz +  Be)  H}|  (2-7) 

Along  an  unloading  path  Ps  then  becomes 

Ps  =  K>‘  .  pe)  -  (Kffl  -  K^)  u*  (l  -  exp[~>J *'p~S'  -])  (2-8) 

where 

u '  =  u  “  u7 

If  the  solution  of  Equation  2-7  produces  a  value  of 
K/  larger  than  «m,  which  is  presumably  the  maximum  bulk  modu¬ 
lus,  it  is  assumed  that  a  phase  change  has  occurred  and  the 
material  is  constrained  to  unload  along  the  linear  path: 

P  =  K  (w*  +  Be)  (2-9) 

where  K  is  P,,/ (um  -  u2). 

The  definition  of  uz  completes  the  specification  of 
Ps.  Experimental  data  (References  2-6,  2-7,  and  2-8)  show  that 
uz  increases  with  increasing  peak  pressure  until  all  of  tin 
cracks  and  voids  are  apparently  closed.  Further  moderate  increase 
in  pressure  does  not  appreciably  increase  the  permanent  set.  At 
very  high  pressures,  u2  may  again  increase  due  to  phase  changes 
in  one  or  more  of  the  constituent  minerals  in  the  rock.  This 
behavior,  depicted  schematically  in  Figure  2-3,  is  expressed 
mathematical ly  by 
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(1-10) 


ivi  -  \> 

U  +  (P  -  P  )  d  P  <  P  •  -EE—; — -  +  P3 
p  in  a  am  d  a 


y  is  the  excess  compression  at  which  all  the  voids  are  closed. 
The  density  op(  =  o0(l  +  up))  must  be  less  than  or  equal  to  the 
crystal  density  of  the  rock.  Pa  is  the  minimum  pressure  at  which 
a  phase  change  can  occur  and  (  =  pQ(l  +  upp))  is  the  crystal 
density  of  the  rock  after  the  phase  change  is  completed.  The 
assumption  that  u2  varies  linearly  with  pressure  (d  is  the 
proportionality  constant)  is  certainly  crude,  but  the  meager 
unloading  data  and  the  large  scatter  within  that  data  does  not 
justify  further  refinements  at  this  time. 

For  the  highly  porous  type  of  material  whose  loading 
hydrostat  is  illustrated  in  Figure  2-2,  the  pressure  during  unload¬ 
ing  is  found  by  evaluating  Equations  2-7  and  2-8  wi  th 
u1  «  u  -  u0.  The  result  of  this  procedure  is  an  unloading  hydro¬ 
stat  which  passes  through  the  point  (Pm,  ym)  and  is  parallel  to 
the  master  curve.  Hysteresis  is  automatically  accounted  for 
during  the  loading  phase. 

In  using  this  model  in  practical  calculations,  mean 
tensile  stresses  may  develop.  Because  of  the  low  value  of  tensile 
strength  for  in  situ  rock,  this  is  prevented  and  Pg  is  set  equal 
to  zero  for  u  <  uz. 


CALCULATION  OF  0EVIAT0RIC  STRESS 

A  complete  description  of  the  material  in  the  solid  state 
requires  the  specification  of  deviatoric  stresses  (oj  .)  as  well  as 
the  mean  stress,  ’f  the  material  behaves  elastically,  the  stress 
increment  is  calculated  from 


=  2G(dr. ! . ) 


(2-11) 
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where 

de! .  =  Deviatoric  strain  increment 

ij 

G  =  Shear  modulus 


Stress  states  are  calculated  incrementally  to  allow  for  the  possi¬ 
bility  of  plastic  deformations.  A  trial  deviatoric  stress  tensor 
is  calculated  at  the  end  of  each  time  step  from 


^°i trial 


=  (a! .)  ,  ,  +  (do! .) 

ij  old  ij 


(2-12) 


where  daj :  is  obtained  from  Equation  2-11.  A  yield  criterion  is 
used  to  determine  whether  the  devjatoric  stress  (o{j) trial  defines 
an  elastic  state.  If  so,  then  (ojj)trial  's  the  correct  devia¬ 
toric  stress,  otherwise,  plastic  deformation  has  occurred  during 
the  time  step,  and  do;:  must  be  adjusted  according  to  the  yield 
criterion  and  the  associated  plastic  potential  flow  rule. 


The  mathematical  model  of  the  shear  modulus  G  is  simi¬ 
lar  to  that  of  the  bulk  modulus.  G  varies  from  an  initial  value, 
G0,  at  normal  density,  to  the  intrinsic  value,  Gmax,  as  the 
pores  and  cracks  are  closed. 


G  -  G  -  (G  -  G  )  exp  ~ 
max  max  o  r  p- 


(2-13) 


Some  experimental  data  (References  2-9  and  2-10)  indicate  that,  for 
rock  which  has  been  cracked  prior  to  testing,  GQ  is  close  to 
zero.  (Reference  2-6  reports  GQ  -  0  ±  k  kb  for  fractured  UTS 
granite.)  However,  as  u  increases,  G  approaches  its  intrinsic 
value  much  more  rapidly  than  does  K.  We  assume  that  variability 
of  the  effective  bulk  modulus  is  strongly  affected  by  the  closing 
of  spherical  pores,  whereas  the  effective  shear  modulus  is  more 
influenced  by  the  closing  of  in  situ  cracks.  It  is  consistent  with 
the  findings  of  Reference  2-3,  which  shows  that  cracks  are  more 
easily  closed  than  pores,  for  G  to  approach  its  intrinsic  value 
faster  than  does  K. 

Specifying  that  G  increases  with  increasing  p  raises 
the  possibility  that  energy  might  be  extracted  from  the  material 
by  hydrostatically  compressing  it,  shearing  it  at  high  pressure, 
releasing  the  pressure,  and  then  releasing  the  shear.  This  danger 
is  avoided  by  assuming  that  friction  prevents  cracks  from  reopening 
during  unloading  so  that  the  largest  value  of  G  reached  on  load¬ 
ing  is  retained  during  subsequent  unloading/ reloading.  Under  these 
restrictions,  a  material  may  dissipate  energy  in  shear  during  load¬ 
ing  and  unloading  cycles  but  can  never  produce  additional  energy. 
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Data  on  the  temperature  dependence  of  material  strength 
(Reference  2-11)  clearly  indicate  that  the  shear  strength  of 
several  rocks  decreases  almost  linearly  with  increasing  tempera¬ 
ture.  By  analogy,  a  similar  dependence  is  assumed  for  the  shear 
modulus.  As  the  specific  internal  energy  approaches  e  , 

G  tends  towards  zero  according  to  the  following 
equation : 


G  = 


max 


(G 


max 


G  )  exp 
o  r 


ill1  - 1) 


(2-14) 


For  e  >  e_,,  the  rock  is  assumed  to  have  undergone  a  phase  change 
such  as  melting,  and  to  be  unable  to  support  shear  stress,  in 
which  case  the  shear  modulus,  G,  is  set  to  zero. 

The  yield  criterion  determines  the  maximum  deviatoric 
stress  which  the  material  can  support  before  it  deforms  plastically. 
The  yield  criterion  is  a  function  of  stress  components.  If  the 
criterion  is  not  satisfied,  the  material  is  assumed  to  behave 
elastically.  If  the  trial  deviator  stresses  calculated  by  Equa¬ 
tion  2-11  are  in  the  forbidden  region  outside  the  yield  surface, 
the  stress  state  is  adjusted  by  means  of  the  flow  rule  so  as  to 
be  exactly  on  the  yield  surface.  The  yield  surface  used  is  a  com¬ 
bination  of  the  llohr-Coulomb  and  von  Mises  representations 
(Reference  2-12),  i.e., 


y[^2 


=  Y  =  minimum  of 


k,  +  k2P 


(2-15) 


where 


J2 


ki 

k2 

k3 


Second  invariant  of  the  stress  deviator;  a 
function  of  the  current  state  of  stress 

Cohesion 

Tangent  of  the  angle  of  internal  friction 
Upper  limit  of  yield  strength  (von  Mises  Surface) 


Towles  and  Riecker  (Reference  2-11)  have  shown  that  the 
shear  strength  of  several  rocks  decreases  with  increasing  tempera¬ 
ture  roughly  in  accord  with  the  following  empirical  relation: 


Y  =  Y  exp  (tT/T  ) 
o  m 


(2-16) 
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where  t  is  a  negative  empirical  constant  whose  value  for  the 
rocks  tested  lies  between  -0.87  and  -1.4,  and 

T  =  Temperature 

Tm  =  Melting  temperature 

Yq  =  Shear  strength  at  room  temperature 

To  include  temperature  effects  in  the  model,  the  exponential  factor 
in  Equation  2-16  is  represented  by  the  first  two  terms  of  its 
Taylor  series  expansion,  t  is  assumed  to  equal  -1,  and  T/Tm  is 
replaced  by  e/e  ,  assuming  constant  heat  capacity.  Thus, 

Equation  2-15  becomes 

__  (k1  +  k2P 

*/j'  =  Y  =  minumum  of  <  (2-17) 

w  ■  y 

As  the  specific  internal  energy  e  approaches  the  melting  energy 
em  and  the  solid  approaches  the  fluid  state,  the  stress  deviators 
are  reduced  to  zero  by  both  the  shrinking  of  the  von  Mises  portion 
of  the  yield  sur.ace  and  the  reduction  of  the  shear  modulus  G. 

The  plastic  potential  flow  rule,  used  to  calculate  devia- 
toric  stress  increments  from  strain  increments  when  the  material  is 
deforming  plastically,  can  be  summarized  by  the  statement  that  among 
all  stress  states  which  lie  on  the  yield  surface,  the  stress  actually 
reached  in  a  given  plastic  strain  increment  is  one  for  which  the 
plastic  work  increment  is  stationary.  The  mathematical  formulation 
of  this  flow  rule  is  discussed  more  fully  in  Reference  2-13.  This 
flow  rule  requires  the  separation  of  the  strain  tensor  into  elastic 
and  plastic  components  and  can  result  in  permanent  plastic  volume 
change.  To  calculate  Ps  after  plastic  volume  change  has  occurred, 
the  elastic  component  of  excess  compression  u  is  used  instead  of 
the  total  excess  compression,  where 

- ! —  =  — ! - -J —  (2-18) 

pe  +  1  u  +  1  u  +  1 

p  ^ 

and  y  is  the  component  of  excess  compression  arising  from  the 
plastic  flow. 
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SECTION  3 

HUG0NI0T  AND  RELEASE  ADIABAT  DATA 


This  section  summarizes  Hugoniot  and  release  adiabat  data 
for  a  number  of  geologic  materials  whose  equation  of  state  has  been 
formulated  as  described  above.  The  materials  considered  and  the 
coefficients  in  the  equations  are  given  in  Table  3"1. 

Comparisons  between  the  equations  of  state  and  the  rele¬ 
vant  data  are  given  in  Figures  3-1  through  3“21 .  For  some 
materials,  comparison  is  made  only  between  data  and  model  along  the 
principal  Hugoniot,  while  for  others,  data  and  models  for  release 
adiabats  are  also  included.  For  these  materials  where  release 
adiabat  data  are  lacking,  the  model  release  adiabat  is  shown  cen¬ 
tered  at  200  kb,  400  kb,  600  kb,  and  800  kb. 
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PRSCEC  fITS  * 


METHOD  OF  FITTING  THE  PRESENT  MODEL 
TO  HUGONIOT  AND  RELEASE  ADIABAT  DATA 


The  user  of  the  equations  of  state  presented  above  may 
wish  to  alter  the  recommended  coefficients  in  order  to  represent 
similar  materials  under  different  conditions  of  initial  density  or 
moisture  content.  He  may  also  wish  to  derive  coefficients  directly 
from  data  for  entirely  different  materials.  The  following  dis¬ 
cussion  of  the  model  for  NTS  granite  is  intended  to  help  the  new 
user  do  this  by  indicating  how  various  material  parameters  are 
evaluated.  Following  the  discussion  of  fitting  a  complete  model 
for  NTS  granite,  some  comments  on  fitting  a  model  to  release 
adiabat  data  are  made. 


DERIVATION  OF  A  MODEL  FOR  NTS  GRANITE 

The  density  of  granitic  rocks  may  vary  from  2.5  to 
2.8  gm/cm^,  while  samples  from  the  region  of  interest.  Area  15  of 
the  Nevada  Test  Site  (NTS),  vary  between  2.62  and  2.70  gm/cm^. 
Densities  are  usually  determined  from  measurements  on  small, 
competent  samples  and,  hence  represent  an  upper  limit  to  the 
average  density  of  the  rock  in  a  volume  comparable  to  that  in  the 
much  larger  finite  difference  zones.  With  this  in  mind,  an 
initial  density  of  2.65  gm/cm^  was  selected  to  represent  the 
granodiorite  under  consideration. 

The  calculation  of  pressure  Pg  for  granite,  a  solid 
exhibiting  hysteresis,  requires  the  evaluation  of  nine  constants 

8,  Kmax»  Ko»  u*»  c’  d’  V  UPP’  atld  Pa‘  The  f'rst» 
is  the  coefficient  of  volume  expansion  per  unit  internal  energy, 

and  can  be  derived  from  the  ratio  of  the  coefficient  of  volume 
expansion  per  degree  Centigrade,  and  the  specific  heat.  Using 
values  of  25  pin. /in.  and  0.192  cal/gm  (Reference  A-l)  for  the  two 
above  mentioned  quantities,  g  =  3.0  (cm/ cm)  (Mb-cc/gm)  . 

Although  neither  the  specific  heat  nor  the  coefficient  of  volume 
expansion  is  constant  over  the  range  of  pressure  and  temperature 
of  interest,  both  tend  to  increase  linearly  with  temperature,  and 
it  is  assumed  here  that  their  ratio  does  not  change.  A  more  elab¬ 
orate  representation  of  g  is  unwarranted,  since  ge  in 
Equation  2-5  is  usually  small  relative  to  the  compression  unless 
the  material  is  subjected  to  pressures  on  the  order  of  hundreds  of 
ki lobars. 
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The  values  of  Kmax,  Kg,  and  J:  for  UTS  granite  can 
be  determined  directly  from  experimental  data  for  pressures  below 
100  kb.  Using  a  hydrostatic  press,  Stephens  (References  4-2  and 
4-3)  has  obtained  loading  and  unloading  pressure/volume  data  on 
several  UTS  rocks  for  pressures  up  to  40  kb.  In  addition,  La  Mori's 
measurements  (Reference  4-4)  of  the  hydrostats  of  Westerly  granite 
and  tonal ite,  and  Birch's  data  (Reference  4-5)  on  the  congressional 
wave  speed  in  various  granites  prestressed  up  to  10  kb,  provide  use¬ 
ful  data  on  possible  variations  in  the  hydrostat.  Triaxial  com¬ 
pression  tests,  in  which  all  of  the  stress/strain  components  have 
been  measured  (References  4-6  and  4-7)  ,  have  been  performed  at 
pressures  less  than  1  kb,  and  these  provide  a  further  check.  In 
an  attempt  to  correlate  these  experiments,  aP/du  was  calculated 
from  the  data  reported,  and  plotted  against  P.  Curves  represent¬ 
ing  the  data  obtained  are  compared  with  the  model  in  Figure  4-1. 

The  variations  among  the  measurements  are  sc  large  that  it  is  not 
possible  to  choose  a  single  value  for  the  "intrinsic"  bulk  modu¬ 
lus,  Kmax,  for  all  granites.  While  all  of  the  curves  tend 
toward  values  between  0.6  and  1.0  Mb,  the  data  vary  too  much  to  be 
fit  by  a  single  value  of  Kmax-  Differences  in  mineral  content 
are  probably  responsible  for  the  measured  variations  in  this 
parameter.  A  value  of  0.8  Mb  is  chosen  as  representative  of  the 
most  recent  data  on  UTS  granodiorite  (Reference  4-3).  Similarly, 
although  KQ  ranges  between  0  and  0.7  Mb  for  the  variety  of 
granites  examined,  a  value  of  0.225  Mb  appears  appropriate  for  the 
UTS  material . 


The  parameter  u*  determines  the  rate  at  which  the  bulk 
modulus  rises  from  KQ  to  l^nax>  while  the  constant  c  is  used 
to  decrease  the  bulk  modulus  at  high  temperatures  to  match  the 
Hugoniot  data  in  the  transition  region  between  solid  and  fluid. 
Values  of  0.0375  and  0.35  for  these  two  parameters  complete  the 
description  of  the  loading  hydrostat.  In  Figure  4-2,  calculated 
hydrostat  and  Hugoniot  curves  are  compared  with  data  for  several 
granites  for  pressures  in  the  range  0  to  45  kb.  Since  the 
Hugoniot  measurements  are  performed  on  small  samples  of  highly 
competent  material,  most  of  those  data  lie  above  the  curve  calcu¬ 
lated  for  the  in  situ  material.  However,  the  upper  curve,  which 
was  calculated  by  assuming  a  constant  shear  modulus  appropriate  to 
consolidated  granite,  fits  the  data  well  up  to  35  kb.  At  higher 
pressures,  the  measurements  on  Westerly  granite  fall  below  the 
calculation.  This  is  not  unexpected  as  the  bulk  modulus,  (curve 
4a,  Figure  4-1),  for  Westerly  granite  and,  therefore,  its  hydro¬ 
stat,  are  also  less  than  those  of  NTS  granite  at  these  pressures. 

The  parameters,  d,  i.  ,  u  ,  and  Pa  are  used  in 
describing  the  unloading  behavior  of  solid  material.  Walsh 
(Reference  4-8)  has  shov/n  that  the  pressure  necessary  to  close  a 
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CURVE  HATER  I AL  REF 

1  MTS  GRANODIGRITE  U-20 

2  NTS  GRAMODIORITE  1,-3 

3a, b  CEDAR  CITY  TONAUTE  . 

Aa  WESTERLY  GRANITE  2A°C  J 

Ab  WESTERLY  GRANITE  160°C  1*'* 

Ac  WESTERLY  GRANITE  3I3°C  ) 

5  ROCKPORT  GRANITE  ) 

6  STONE  MT.  G  WESTERLY  GRANITE/ A-5 

7  PORTERVILLE  GRANITE  I 

8  STONE  HT.  GRANITE  1 

9  WESTERLY  GRANITE  /A-7 

•  NTS  GRANITE  MODEL  7 


O.OI6  0.020  0.02A  0.028  0.032  0.036 


PRESSURE,  MB 


AJAOIA 


PRESSURE  (LOADING)  FOR  GRANITIC  ROCKS 


A 
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crack  is  linearly  related  to  the  dimensions  of  the  crack.  If 
there  are  many  randomly-sized  small  cracks  which  contribute  to  the 
initial  porosity,  it  is  reasonable  to  expect  an  irreversible  change 
in  porosity  which  varies  inversly  with  peak  loading  pressure  until 
all  cracks  are  closed.  Mo  further  permanent  change  in  porosity  is 
then  exhibited  until  the  pressure  is  so  high  that  either  the  pores 
break  down  or  an  irreversible  phase  change  occurs.  This  behavior 
is  modeled  by  Equation  2-10,  where  d  is  the  constant  of  propor¬ 
tionality  between  the  loss  of  porosit/  and  the  pressure;  p  is 
the  total  change  in  compression  when  all  the  cracks  are  cioied; 

Pa  is  the  pressure  above  which  pore  breakdown  and/or  phase  change 
is  possible;  and  upp  is  the  maximum  irreversible  change  in  com¬ 
pression  due  to  pore^breakdown  or  phase  change.  Unloading  paths 
from  the  granite  Hugoniot  are  shown  in  Figure  3~2.  At  pressures 
above  400  kb,  enough  internal  energy  is  present  in  the  material 
so  that  the  fluid  pressure  P,,  is  not  zero.  The  unloading  path 
from  these  high  pressures  drops  very  sharply  until  the  contribution 
of  Pf  to  the  total  pressure  dominates.  Beyond  that  point,  the 
pressure  diminishes  so  slowly  that  the  material  may  not  actually 
return  to  zero  pressure  until  it  has  expanded  beyond  the  initial 
sol i d  dens i ty . 

Evaluating  the  parameters  for  Pf  in  Equation  2-2  is 
all  that  remains  to  define  the  behavior  of  the  mean  stress  in  MTS 
granite.  At  low  densities,  the  exponential  expression,  exp  2, 
damps  out  everything  but  the  ajpe"  term.  Hence,  aj  is  the 
equivalent  of  y  -  1  in  the  perfect  gas  law.  Calculations  at 
Lawrence  Livermore  Laboratory  (Reference  4-9)  show/  that  y  -  1 
asymptotically  approaches  a  value  between  0.04  and  0.14  for  several 
rock  media  at  low  densities  for  pressures  below  10  kb.  Consequently, 
aj  was  set  to  0.1  for  all  of  the  rocks  modeled  in  this  study.  The 
coefficient  a  of  u/n  in  the  exponential  was  set  to  5,  consistent 
with  the  studies  on  metals  by  Tillotson  (Reference  4-10)  and  on 
rocks  by  Allen  (Reference  4-11).  For  compressed  material  at 
energies  considerably  greater  than  the  vaporization  energy,  the 
term 


becomes  negligible  and  the  coefficient  of  pe"  is  simply 

aj  +  a 2»  The  term  a2  was  set  at  0.4  so  that  the  sum  equaled  0.5, 

again  consistent  with  Tillotson  and  Allen.  This  model  does  differ 


from  Ti 1 lotson's  formulation  in  that  the  term 


is  set  to  zero 


at  energy  densities  less  than  that  required  to  melt  the  material. 
For  granite,  which  melts  at  a  temperature  of  about  1400°C,  the  mini¬ 
mum  energy  density  for  melting,  e^,  is  assumed  to  be 
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0.0115  Mb-cc/gm.  Since  the  melting  temperature  in  most  materials 
increases  with  compression,  the  energy  density  at  melting,  em 
is  allowed  to  increase  linearly  to  a  maximum  emm,  equal  to 
0.035  Mb-cc/gm  in  granite,  so  that 


e  <  e 
mo  m 


e  (1  -  f n) 
mo 


<  e 

mm 


(4-1) 


where  f  is  set  to  1.7.  In  Figure  3-1  the  calculated  principal 
Hugoniot  is  compared  with  experimental  data  up  to  1  Mb.  Data  for 
several  rocks  with  mineral  contents  similar  to  NTS  granite  are 
included.  Above  100  kb,  all  of  these  rocks  behave  alike,  suggest¬ 
ing  that  their  basic  chemical  composition,  rather  than  initial 
density,  porosity  or  physical  structure,  largely  determines  their 
behavior  at  high  pressures. 

The  mean  stress  behavior  of  NTS  granite  is  now  determined 
and  attention  is  shifted  to  the  parameters  influencing  its  devia- 
toric  characteristics.  Except  for  Simmons  work  (Reference  4-2)  on 
the  shear  wave  velocities  of  rocks  under  hydrostatic  pressures  up 
to  10  kb,  there  is  little  data  from  which  to  determine  the  shear 
modulus,  G,  of  granite.  Recently  improved  techniques  in  triaxial 
testing  (Reference  4-8)  show  considerable  promise  as  a  means  of 
determining  G,  but  the  work  is  still  preliminary.  Stephens 
(Reference  4-13)  at  Lawrence  Livermore  Laboratory  measured  the 
shear  modulus  for  consolidated  and  cracked  NTS  granodiori te.  His 
data,  the  wave  speed  measurements  of  Simmons,  and  the  model  are  pre¬ 
sented  in  Figure  4-2.  The  values  of  GQ,  Gmax,  and  y*  match 
Stephens'  data  for  the  cracked  granodior i te. 


Data  on  the  yield  strength  of  granite  as  a  function  of 
mean  stress  are  presented  in  Figure  4-3.  Several  Mohr-Coulomb 
surfaces  are  also  plotted  for  comparison.  The  large  differences 
between  the  various  experiments  are  attributable  to  the  initial 
condition  of  the  rock,  since  preexisting  cracks  and  pores  or  varia¬ 
tions  in  the  water  content  and  pore  pressure  can  severely  affect  the 
strength  of  the  rock.  There  is  evidence  that  the  yield  surface 
depends  to  some  extent  on  the  stress  state  (References  4-14  and 
4-15),  but  this  is  not  included  in  this  model.  For  intact,  compe¬ 
tent  granite,  the  data  is  best  fit  by  choosing  kj  =  0.0005  Mb 
and  k£  =  1.10.  The  measurements  made  on  jointed  and  cracked 
samples  from  the  Pile  Driver  site  in  NTS  Area  15  lead  to 
kj  =  0.0001  Mb  and  0.5  <  k£  <  0.7  (References  4-6  and  4-16). 
Finally,  recent  unpublished  Lawrence  Livermore  Laboratory  data 
(Reference  4-17)  from  presumably  cracked,  wet  NTS  granite,  suggest 
that  k£  may  be  as  low  as  0.3. 
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Since  the  small  samples  used  in  the  experiments  are 
generally  more  competent  than  in  situ  material,  the  measured  yield 
strengths  have  been  treated  as  upper  bounds  to  the  values  to  be 
expected  in  the  field  and  values  of  kj  =  0.0001  Mb  and 
k2  =  0.3  are  used  in  the  model.  For  similar  reasons,  the  von  Mises 
surface,  k^,  is  set  at  0.017  Mb,  although  comparison  of  the  data 
near  the  Hugonioi  elastic  limit  with  the  calculated  hydrostat 
suggests  a  value  of  as  high  as  0.020  Mb.  The  model  at  low 

pressure  consists  chiefly  of  the  hydrostat,  shear  modulus,  yield 
criterion  and  flow  rule.  These  are  checked  by  comparing  the  model 
hydrostat  and  hugoniot  with  the  data  in  Figure  4-4. 


APPLICATION  OF  THE  MODEL  TO  MATCH  RELEASE  ADIABAT  DATA 

Fitting  the  model  to  release  adiabat  data  is  performed 
by  selecting  the  parameters  u,  Pg,  upp,  e  ,  and  f. 

The  amount  of  hysteresis  which  occurs  before  melting, 
that  is  when 


e  <  em 
-  m 


is  controlled  through  the  parameters  u  ,  up  ,  Pa,  and  d  as 
indicated  in  Equation  2-10.  Increasing^  up  P£nd  ppp  increases 
the  amount  of  hysteresis  and  the  slope  of  the  release  adiabat. 


Behavior  such  as  that  shown  in  Figure  4-5  is  controlled 
through  the  parameters  e^  and  f.  By  raising  or  f,  the 

energy  required  to  vaporize  the  solid  is  raised  and  the  solid 
material  remains  a  solid  (Case  2).  By  lowering  emo  or  f,  the 
energy  required  to  vaporize  the  solid  is  decreased  and  the  pressure 
may  still  be  governed  appreciably  by  the  gas  contribution;  in  this 
case  when  the  pressure  gets  low  enough,  expansion  beyond  the 
initial  volume  can  occur  (Case  1). 


MODEL  HYDROSTAT 


R-71 3^-2283 

REFERENCES 

4-1.  Handbook  of  Chemistry  and  Physics ,  Chemical  Rubber  Publish¬ 
ing  Company,  36th  Ed.,  1 955 >  pp.  2061-2103. 

4-2.  Stephens,  0.  R. ,  and  E.  M.  Lilley,  Static  PV  Curves  of 
Cracked  and  Consolidated  Earth  Materials  to  40  Kilobars, 
UCRL-14711,  Lawrence  Radiation  Laboratory,  March  1 966 . 

4-3-  Private  communication  from  D.  R.  Stephens  (Lawrence  Radia¬ 
tion  Laboratory),  1 968. 

4-4.  LaMori,  P.  N. ,  Compressibility  of  Three  Rocks,  (A)  Westerly 
Granite  and  Solenhofen  Limestone  to  40  Kilobars  at  300°C, 

(B)  Cedar  City  Tonalite  to  40  Kilobars  at  Room  Temperature, 
DASA-2151 ,  Defense  Atomic  Support  Agency,  August  1 968 . 

4-5.  Birch,  F.,  "Velocity  of  Compress ional  Waves  iii  Rocks  to 

70  Kilobars,  Part  1,"  J.  Geophys.  Res.,  Vol.  65*  No.  4,  i960, 
p.  1083. 

4-6.  Tests  for  Strength  Characteristics  of  Rock,  Pile  Driver 

Project,  MRDL  64/90,  U.  S.  Army  Corps  of  Engineers,  Missouri 
River  Division  Laboratory,  1964. 

4-7.  Brace,  W.  F.,  "Some  New  Measurements  of  Linear  Compressibil¬ 
ity  in  Rocks,"  J.  Geophys.  Res.,  Vol.  70,  No.  2,  1965,  p.  391. 

4-8.  Walsh,  J.  B.,  "The  Effect  of  Cracks  on  the  Compressibility 

of  Rock,"  J.  Geophys.  Res.,  Vol.  70,  No.  2,  January  15,  1965. 

4-9.  Butkovitch,  T.  R. ,  The  Gas  Equation  of  State  for  Natural 
Materials,  UCRL-14729,  Lawrence  Radiation  Laboratory, 

January  1 967. 

4-10.  Tillotson,  J.,  Metallic  Equations  of  State  for  Hypervelocity 
Impact,  GA-3216,  General  Atomic  Div.,  General  Dynamics 
Corporation,  July  1 962 . 

4-11.  Allen,  R.  T. ,  Equation  of  State  of  Rocks  and  Minerals, 

GAMD-7834,  General  Atomic  Div.,  General  Dynamics  Corpora¬ 
tion,  1967. 

4-12.  Simmons,  G.,  "Velocity  of  Shear  Waves  in  Rocks  to  10  Kilo¬ 
bars,  Part  1,"  J.  Geophys.  Res.,  Vol.  69,  No.  6,  1964,  p.  1123. 

4-13.  Stephens,  D.  R. ,  Elastic  Constants  of  Fractured  Granodiorite , 
UCID-I5369,  Lawrence  Radiation  Laboratory  (to  be  published). 


43 


vnMVMnMv**  „ 


4-14.  Mogi,  K. ,  "Effect  of  the  Inu  mediate  Principal  Stress  on 
Rock  Failure,"  J.  Geophys.  Res.,  Vol.  72,  No.  20,  October 
1967,  p.  5P7. 


t 

* 


i 


1 

t 

t 

i 

i 

1 


4-15.  Brown,  W.  S.,  et  al.,  Influence  of  Dynamic  Loading,  Biaxial 
Loading  and  Prefracturing  on  the  Stress-Strain  and  Fracture 
Characteristics  of  Rocks,  DASA-2713,  University  of  Utah, 
March  1971 . 

4-16.  Byerlee,  J.  D.  ,  "Frictional  Characteristics  of  Granite 

Under  High  Confining  Pressure,"  J.  Geophys.  Res.,  Vol.  72, 
No.  14,  1967,  p.  3639- 

4-17.  Private  communication  from  H.  Heard  (Lawrence  Radiation 

Laboratory)  via  D.  R.  Stephens  (Lawrence  Radiation  Labora¬ 
tory)  ,  1968. 

4-18.  Brace,  W.  F. ,  et  al.,  "Di  latency  in  the  Fracture  o'  Crystal¬ 
line  Rocks,"  J.  Geophys.  Res.,  Vol.  71,  No.  !6.  August  15, 
1966. 

4-19.  Giardini,  A.  A.,  et  al . ,  "Triaxial  Compression  Data  on 

Nuclear  Explosion  Shocked,  Mechanically  Shocked  and  Normal 
Granodiorite  from  the  Nevada  Test  Site,"  J.  Geophys.  Res., 
Vol.  73,  No.  4,  February  15,  I960,  p.  305. 

4-20.  Stephens,  D.  R.,  and  E.  M.  Li  1  ley,  Static  PV  Curves  of 
Cracked  and  Consolidated  Earth  Materials  to  40  Kilobars, 
UCRL-14711,  Lawrence  Radiation  Laboratory. 

4-21.  Jones,  A.  H.,  and  N.  H.  Froula,  Uniaxial  Strain  Behavior 
of  Four  Geological  Materials  to  50  Kilobars,  DASA-2209, 

March  1 969  - 

4-22.  Grine,  D.  R.,  Equations  of  State  of  Granite  and  Salt, 
UCRL-13004,  Lawrence  Radiation  Laboratory,  May  1961. 

4-23.  Ainsworth,  D.  L.,  and  B.  R.  Sullivan,  Shock  Response  of  Rock 
at  Pressures  Below  30  Kilobars ,  WE S  TR  6-802,  U.  S.  Army 
Engineer  Waterways  Experiment  Station,  November  1 967 • 


44 


f: 

>3 


- 


ak£U.i  -JuiiVi LaS. if. 'ZL „ 


.  • 


R-71 34-2283 


APPENDIX  A 


During  the  course  of  fitting  and  checking  the  models 
described  above,  references  containing  release  adiabat  data  were 
found.  These  references,  some  of  which  themselves  contain  exten¬ 
sive  lists  of  references,  are  given  below. 


TABLE  A-l. 


Pressure  Range, 

Material  k  Bar  Reference 


Alluvium,  Frenchman's  Flat,  NTS 

1 

to 

12 

A- 9 

Playa,  Area  5,  NTS,  0  to  18.92 

Moisture 

70 

to 

280 

A- 2 

Tuff,  Area  12,  NTS,  dry 

85 

to 

125 

A- 4 

Tuff,  Area  19,  NTS,  dry 

20 

to 

290 

A- 8 

Tuff,  Area  12,  NTS,  saturated 

115 

to 

145 

A- 4 

Tuff,  Rainier  Mesa,  NTS,  saturated 

135 

to 

470 

A- 6 

Sandstone,  Coconino,  dry 

155 

to 

255 

A- 4 

Quartz,  Arkansas  Novaculite 

70 

to 

150 

A-l 

Quartz,  Arkansas  Novaculite 

70 

to 

380 

A- 3 

Quartz,  Fused  Glass 

100 

to 

380 

A-6 

Granite,  Raymond,  California 

45 

to 

285 

A- 5 

Granodiori te,  Climax  Stock,  NTS 

190 

to 

285 

A- 8 

Tonal ite.  Cedar  City,  Utah 

10 

to 

65 

A-9 

Anorthosite,  San  Gabriel  Mtns.,  Calif. 

55 

to 

422 

A- 3 

A1  luvium 

80 

to 

285 

A- 11 

45 


k. 
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